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Abstract 

We present a detailed study of the static enhancement effects of electric charges in ^m-long single- 
walled carbon nanotubes, using theoretically an atomic charge-dipole model and experimentally 
electrostatic force microscopy. We demonstrate that nanotubes exhibit at their ends surprisingly 
weak charge enhancements which decrease with the nanotube length and increase with the nanotube 
radius. A quantitative agreement is obtained between theory and experiments. 
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Understanding of the properties of electric charges in carbon nanotubes (CNTs) is one of 
the important issues for their promising apphcations in nanoelectromechanical systems,^ field 
emission^ chemical sensor^ and charge storage .'^'^ A key-aspect of the electrostatics of these 
one-dimensional systems is the knowledge of the distribution of electric charges along the 
nanotubes, because charges are likely to accumulate at the nanotube ends due to Coulomb 
repulsion. Theoretical predictions have been established for this effect, but not in the range 
of lengths accessible from experiments, so that no comparison has been established between 
theory and experimental observations so far. More precisely, the electrostatic properties 
of single-walled nanotubes (SWCNTs) have been addressed on the one hand using electric 
force microscopy (EFM) experiments^ coupled to charge injection technique d^*^" *. Results 
obtained for /im-long nanotubes indicated that electric charges are distributed rather uni- 
formly along the tube length, with however no theoretical support in this range of nanotube 
length. On the other hand, in theoretical studies, density functional theorj^^I and classi- 
cal electrostatic^ calculations have been performed to compute the charge distribution in 
SWCNTs, and have predicted ?7-like shapes due to a charge accumulation at the nanotube 
ends. These calculations however only hold for short (< 100 nm) nanotube lengths, which 
are not easily accessible from experiments. 

It is the scope of this paper to provide a combined experimental and theory work on 
this issue. We present a detailed study of the static enhancement effects of electric charges 
in SWCNTs, using theoretically an atomic charge-dipole model and experimentally electro- 
static force microscopy. It is demonstrated that the U-like shape of the charge distribution 
expected for short nanotubes is replaced in the case of /xm-long tubes by weak charge en- 
hancements localized at the nanotube ends, in agreement with the experimental values for 
the enhancement factors (up to few tens of %) observed from EFM and charge injection 
experiments. The dependence of the charge enhancement factors on the nanotube radius 
has also been measured from EFM experiments, and falls in quantitative agreement with 
theoretical preditions for yum-long tubes. 

The paper is organized as follows : we first describe the numerical calculations of the 
charge distribution along nanotubes using the atomic charge-dipole model and the results 
obtained for short nanotubes with open or closed caps, either considered in vacuum or on 
a Si02 substrate. The extrapolation procedure to the case of /xm-long nanotubes is then 
presented, and compared with experimental EFM measurements of charge enhancement 
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factors on SWCNTs. We finally discuss the dependence of the enhancement factors as a 
function of the nanotube radius. 

In the theoretical calculations presented throughout this work, the interactions between 
the electric charges and the induced dipoles are described using the Gaussian-regularized 
atomic charge-dipole interaction model,'^^'^ in which the atoms are treated as interacting 
polarizable points with free charges, and the distribution of charges and dipoles are deter- 
mined by the fact that their static equilibrium state should correspond to the minimum 
value of the total molecular electrostatic energy. Compared with classical Coulomb-law- 
based models in which only charges are considered, this model provides a more accurate 
description of electrostatic properties of CNTs, since the charges, the induced dipoles and 
the atomic polarizabilities are taken into account. 

In order to achieve a valid comparison between experimental data and calculation results, 
the effect of a Si02 substrate (nanotubes are usually deposited on a Si02 thin film in experi- 
ments) is also taken into account in our calculations, by adding surface-induced terms to the 
vacuum electrostatic interaction tensors using the method of mirror images.'^ The dielectric 
constant of Si02 is taken as 4.0. The average distance between the bottom of the tubes and 
the Si02 surface is set to d = 0.34 nm after the computed CNT-Si02 long-range interacting 
configurations from Refs.'^^'^. Furthermore, we note that d can slightly vary with the tube 
radius R. It will however be fixed to 0.34 nm in this work as an average value. The atomic 
structure of CNTs is then optimized by energy minimization using the method of conjugated 
gradient based on a many-body chemical potential model AIREBO (adaptive interatomic 
reactive empirical bond order) Linear charge densities of 0.055e/nm have been used in 
calculations, so as to match linear charge densities observed experimentally.'^ 

To illustrate the typical outputs of the atomic-scale calculations, we show in Fig.[T] the 
charge distribution at the end of a (9, 0) CNT of length L ^ 11.5 nm and average charge 
density cr'*^'^ = —6.6 x 10~^e/atom. The color of the atoms is proportional to their charge in 
the figure. We represented here for sake of clarity the four distinct situations in which the 
nanotube exhibits either a closed (Fig.jl^ and c) or an open (Fig.Q) and d) cap structure, 
and the tube is either considered in vacuum (Fig.jl^ and b) or deposited on a Si02 substrate 
(Fig.jlj:; and d). As seen from Fig.[l| the charge density at the tube ends is higher than that 
at other parts of the tubes in all situations. The maximum charge density on the opened 
cap is here about twice that on the closed one for this small-radius tube. Finally, when the 
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FIG. 1: (Color online) Charge distribution at the ends of an open-ended and a closed-ended (9, 0) 
SWCNTs (L w 11.5 nm, cr™'= = -6.6 x 10~^e/atom) in free space ((a) and (b)) and on the Si02 
substrate((c) and (d)). The color of the atoms is proportional to the local charge density. The green 
vectors stand for the induced atomic dipoles. The dark arrows stand for the local electric fields 
induced by the net charge, their length and color are proportional to the field intensity, (a) The 
minimum and maximum atomic charge densities in this tube are: o"™" = —5.3 x 10~^ e/atom and 
^max ^ X -^Q-^ e/atom, respectively, (b) a™*" = -5.2 x 10"^ e/atom and a™^^ = -34 x 10"^ 
e/atom. (c) cj™*" = -h5.6 x 10"^ e/atom and cj™"^' = -34 x 10"^ e/atom. (d) a™" = +5.8 x 10"^ 
e/atom and cr'""^ = —74 x 10~^ e/atom. 



tube is deposited on the Si02 surface, electrons are attracted by their image charge towards 
the Si02 surface, ClS db typical semi-space effect. 

Since the nanotubes used in experiments have lengths in the micrometer range, and since 
this scale can hardly be directly addressed by calculations using atomic models due to the 
limit of computational resources, the issue about the relationship between the tube length L 
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FIG. 2: Charge profile along three (9, 0) SWCNTs with different tube lengths L, in space (hollow 
symbols) and upon a Si02 surface (solid symbols), using a separation distance d = 0.34 nm (see 
text). The total net charge density on each tube o"™*^ is fixed to 6.4 x 10~^e/atom (equivalently, 
0.055e/ nm). Each point is calculated as the average value of the charge carried by the nanotube 
atoms over 10% L. 

and the charge distribution needs to be carefully addressed, so as to later extrapolate charge 
enhancement factors to the length scales of interest in experiments. The length dependence 
of the charge enhancements at the nanotube ends is illustrated in Fig.[2| in which we plotted 
the local average charge density as a function of the position along the nanotube (the x-axis 
origin in Fig. [2] corresponds to the nanotube midpoint). The local average charge density is 
defined from the charge carried by individual CNT atoms, when averaged along the nanotube 
circumference and along a fraction of the length L of the CNT (this fraction is taken as 10%L 
in Fig.[2]). a"''"'^ is the quantity which can be accessed experimentally from EFM techniques.'^ 
The typical shape of the CNT charge distribution observed in Fig. [2] corresponds to the U- 
like shape expected for short nanotubes,'^ but the charge enhancement at the tube ends is 
already seen to become less significant when the tube gets longer. Furthermore, we can also 
see that the charge enhancement is weaker when the ends of the tubes are closed and when 
the nanotubes are placed on the Si02 surface. The latter effect can be understood by the 
fact that the net nanotube charge is located at the CNT side close to the substrate (see 
Fig.[l] (c) and (d)), which leads to an effective reduction of the charge-distributed area in 
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the non-axial direction, similar to an effective decrease of the nanotube radius R, which will 
be discussed further in this paper (see Fig.|6]). 
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FIG. 3: Ratio of charge enhancement 99 as a function of tube length L (in common logarithmic 
scale). This ratio is calculated for both an open-ended (circles) and a closed-ended (triangles) (5,5) 
SWCNTs (radius R 0.34 nm) placed upon a Si02 surface (solid symbols) with d = 0.34 nm, 
and is compared with that for the same tubes in space (empty symbols). The symbols present the 
calculated points, and the lines stand for the extrapolation curves. 

We now focus on charge enhancements for /im-long nanotubes, and their comparison 
with experimental results. Since the spatial resolution in EFM experiments is about lOOnm 
(this resolution is mostly limited by the tip-substrate separation during EFM detection), 
we now consider the enhancement zone in our calculation as a zone of length 10%L at 
the tube end, and define the charge enhancement ratio as the ratio between the charge 
density a'^"''^ averaged in the zone of length 10%L at the end of the nanotube, and the charge 
density a^^'^'^^^ at the center of the nanotube. The influence of the tube length on the charge 
enhancement ratio ip = (j'^^'i / (j"^^d.die -g g];^Q-^n Fig.jsJ ip is seen to decrease significantly 
with L for short tubes (particularly for L < 10 nm), but the variations get smaller when 
the tube is longer. Note that ip is independent of a"'"^, because the local charge densities 
should be proportional to the total one by requiring a constant electric potential on the 
tube surface. Furthermore, we find that if the (open or closed) cap structure plays an 
important role in the charge enhancement for short tubes (L < 15 nm) (as seen in Fig.jl] 
), this effect already becomes unsignificant for L ^ 30 nm, and will become negligible for 

6 



/iin-long nanotubes in experiments with ^ lOOnm resolution. Finally, it appears that the 
only parameter that needs to be properly taken into account is the presence of the Si02 
surface below the nanotube, which still effectively reduces the charge enhancement ration ip 
for L = 30 nm. 

In order to extrapolate these results towards /xm-long nanotubes, we performed a fit 
of the data points of Fig.|3| Since the the analytical formula of the exact distribution 
of charge on a hollow tube is not known in the literature, we used the equation : ip = 
ln(ai X L + a2)/ln(a3 x L + a^), in which a„ (1 < n < 4) are four fitted parameters for each 
nanotube radius. This phenomenological equation has been chosen since it describes a ratio 
between two cylindrical capacitances, and is thus well-suited to account for ip, which is the 
ratio between the linear charge densities at the end and at the middle of the nanotube. The 
lines in Fig.|3] correspond to the fits obtained independently for the nanotubes with either 
open or closed caps, vacuum environment, or Si02 surface. The extrapolated values for ip 
are seen to converge for large L for open and closed cap structures, but to differ depending 
on the vacuum or Si02 environment. This behaviour is in full agreement with the trend 
observed on the atomic calculation points obtained for L 30 nm, which already brings 
confidence at this stage about the validity of our extrapolation procedure. 

Our theoretical predictions are finally compared with electrostatic measurements per- 
formed by injecting and detecting charges in individual CNTs using electrostatic force mi- 
croscopy. In these experiments, nanotubes grown by chemical vapour deposition are de- 
posited from dichloromethane solutions onto silicon wafers covered by a 200nm-thick ther- 
mal dioxide layers. Individual nanotubes are located by atomic force microscopy, and then 
charged (see Fig.|4](a)) by pressing the biased tip of an atomic force microscope on the nan- 
otube (typically with an injection bias Vinj = — 5V, pressing force of a few nN). The CNT 
charge state is then measured before and after injection by EFM, by recording electrostatic 
force gradients acting on the tip which is intentionally lifted at a distance z about 50-100 
nm above the sample surface to discard short range surface forces. Fig.|4] (b) shows the 
topography image of a SWCNT. In Fig.|4] (c), the EFM scan of the tube before charging 
is shown, as a dark footprint of the CNT topography associated with attractive forces due 
to the nanotube capacitance. It can be shown experimentally that the negative frequency 
shifts are here of capacitive origin, and not originating in a positive charge transferred from 
the substrate to the nanotube (see details in Rei^. The nanotube EFM image after charge 
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FIG. 4: (Color online) (a) Schematics of the charge injection and detection with the EFM tip. 
Charging takes place when the biased tip is put in contact with the CNT. During the data acqui- 
sition cantilever is lifted at distance z above the surface, (b) Topography image of a SWCNT with 
0.5 nm radius deposited on 200 nm oxide layer, (c) EFM image acquired before injection. The dark 
feature corresponds to the uncharged tube, (d) EFM scan after injection experiment. The bright 
feature corresponds to the charged tube with the uniform linear charge density of 0.055e/nm. 

injection is shown in Fig.|4] (d). The tube is seen here as a bright feature as a result of 
the negative charges injected in the tube. From previous EFM studies, we have shown that 
the charge imaged for SWCNTs mainly correspond to charge emitted from the tube and 
"printed" in the oxide layer in the vicinity of the nanotube.'^^'^ 

To compare these predictions with our calculation results, we show in Fig.|5]the charge 
distribution at the end of a SWCNT (total length 2 /im) after a charge injection experiment. 
A non-linear color scale has been used in Fig.[5]D in order to evidence the weak charge 
enhancement localized within 200 nm at the nanotube end. The charge distribution along the 
nanotube is shown in Fig.jS}:;, in which we plotted the charge enhancement factor measured 
from EFM, defined as the ratio of the EFM signal with that measured at the middle of the 
nanotube. From these experimental data, one gets the maximum value <y9 = 1.17 ± 0.05 for 
this tube (see Fig.js] (c)), in agreement with the numerical extrapolation from theoretical 
results predicting ^9 ^ 1.165 for a 2 yum tube with R = 0.8 nm deposited on a Si02 surface. 

To further validate the comparison of our theoretical predictions with experiments, we 
now focus on the dependence of the charge enhancement ratio as a function of the nanotube 
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FIG. 5: (Color online) (a) AFM topography image of single nanotube with R = 0.8 nm and L ~ 2 
^m deposited on silicon dioxide, (b) EFM scan of the same tube made after charge injection. A 
non-linear color scale has been used in order to clearly show the weak enhancement at the tube 
end. The black line is a guide-to-the-eye for the physical end of the tube, (c) Experimental charge 
enhancement along the axis of the nanotube, defined as the ratio of the EFM signal with that 
measured at the middle of the nanotube. 

radius R. Such an analysis would not be possible for short nanotubes, because the charge 
enhancement ratio would then be strongly dependent on the nanotube cap structure, as 
discussed previously (see Fig.js]), while this effect is not relevant for /im-long naontubes. 
Intuitively, one can guess that the nanotube charge enhancement factor will increase with 
the tube radius R, because the enhancement factor decreases with the nanotube length L : 
increasing R at fixed length L reduces the nanotube anisotropy, and is qualitatively similar 
as decreasing the nanotube length L for a fixed radius R. 

Experimentally, we measured the charge densities along seven SWCNTs with lengths 
between 1 and 9 /im in a similar way as in Fig.|5] (c), and plotted the corresponding charge 
enhancement ratios in Fig.[6| as a function of the nanotube radius R measured from atomic 
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FIG. 6: Ratio of charge enhancement for a number of tubes with different radius i?, on a 
Si02 surface. The soUd squares stand for derived from the experimental measurements of seven 
nanotubes (with lenght L = 1 ~ 9 /im) deposited on 200 nm sihcon oxide layer. The symbols "+" 
stand for the extrapolation results for ^m length tubes. 

force topography images. The ±0.05 error bars on y9 correspond here to the accuracy of the 
EFM measurements. Experimental data points clearly show that slightly increases as a 
function of the nanotube radius R. The possibility to observe this behavior also confirms 
that the values of on /im-long nanotubes do not critically depend on the tube length, nor 
on the nanotube cap structure. Numerical calculations for the charge enhancement ratio 
(y9 have also been performed using nanotubes with different radius i?, and are shown in 
Fig. [6] Direct calculations of obtained from the atomic dipole-charge models and using 
an averaging over 10% L are given in Fig.[6] for two short nanotubes (9 and 12 nm, solid 
circles and triangles), as well as calculation results obtained for /im-long nanotubes (dotted 
line) using the extrapolation procedure described in Fig.|3} Theoretical predictions are seen 
to quantitatively agree with experimental data within experimental error bars, and confirm 
the increase of the charge enhancement ratio as a function of the nanotube radius. The 
values of computed with 5%L and 15%L for /xm length tubes (data not shown) also vary 
within experimental error bars. 

In summary, we have characterized the enhancement of net electric charge in SWCNTs 
by both atomic-model calculations and EFM experiments. We have demonstrated that 
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the U-like shape of the charge distribution expected for short nanotubes is replaced for 
/xm-long nanotubes by weak charge enhancements locahzed at the nanotube ends, while 
the nanotube charge densities are otherwise almost constant along the nanotubes. The 
dependence on the tube length, nanotube cap structure, and the influence of silica substrate 
have been investigated. It has been shown that the charge enhancement at the ends of 
CNTs depends strongly on the geometry of the cap only for short tubes (< 100 nm), but 
has an insignificant influence for nanotubes with lengths in the micrometer range. The 
increase of the charge enhancement ratio with the nanotube radius has been demonstrated 
experimentally, in quantitative agreement with theoretical predictions. 

We expect that the mapping and the understanding of the charge enhancement of CNTs 
are important for many applications, besides the fundamental character of this study, e.g. 
imaging of field and charge in CNTs electronic circuitd^^ES or emission devices modified 
by the presence of surfaces. The electrostatic response of nanotubes appears to be strongly 
sensitive to its environment, which is of high importance for nanotube based sensors .'^'^ This 
work can also have implication in the field of nano-electromechanical systems and charge 
storage devices. 

This work is a part of the CNRS GDR-E No. 2756. Z. W. acknowledges the support 
from the region of Franche-Comte (grant 060914-10). 
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